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We present a new procedure for improving the effective potential by using 
renormalization group equation (RGE) in the presence of several mass scales. We 
propose a modification of the mass- dependent (MD) renormalization scheme, MD 
scheme, so that the scalar mass parameter runs at most logarithmically on the one 
hand and the decoupling of heavy particles is naturally incorporated in the RGE's 
on the other. Thanks to these properties, the procedure in MD scheme turns out 
to be very simple compared with the regionwise procedure in MS scheme proposed 
previously. The relation with other schemes is also discussed both analytically and 
numerically. 



1. Introduction 



Recently, there have been the renewed interests on how to sum up large log- 
arithms in the effective potential, to investigate the standard model and beyond. 
Basically, large logarithms like ln(M//z), which makes the perturbation expansion 
unreliable, appear when one deals with a system possessing large mass scale M 
compared with the scale \i at which one discuss the physics. In this situation one 
considers resumming the perturbation series by using the renormalization group 
equation! 1 ' When one concerns with the functional form of the effective potential, 
one considers its renormalization-group (RG) improvement. This is well-known 
since the work by Coleman and Weinberg for the massless A0 4 theory, although 
the complete description even for the massive A0 4 theory has been given only 
recently! 3 ' 4 ' 51 

In many realistic applications, one often has to deal with an additional mass 
scale m with the hierarchy \i Cm <C M. In the supersymmetric standard model, 
for instance, one can regard /x, m and M as the weak scale, supersymmetry breaking 
scale and unification scale, respectively. When we discuss such a system, we face the 
problem of multi-mass-scales! 6 ' there appear several types of logarithms, ln(M//x) 
and ln(m//x), while we are able to sum up just a single logarithm by using the 
RGE. 

In Ref. [7], one way to improve the effective potential in the presence of multi- 
mass-scales was described in MS renormalization scheme. The point was to make 

ro -jo] 

use of the decoupling theorem and to divide the energy region (region of field 
space) so that in each region, there remains essentially a single log factor. Although 
there is nothing wrong in principle, such regionwise procedure may be cumbersome 
in practice. So it is desirable to have an alternative way to handle multi-mass-scale 
systems. 

In this paper, we propose a simple modification of the conventional mass- 
dependent (MD) renormalization scheme, which we call modified MD scheme (MD 
scheme), and apply it to improving the effective potential in the presence of sev- 
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eral mass scales. Basically in MD scheme, the RG coefficient functions (f3 and 7 
functions) depend on mass parameters and hence the decoupling of heavy particles 

[13] 

is taken into account in the form of RG runnings. In addition, the proposed MD 
scheme has a property that mass parameters run at most logarithmically while 
keeping the 'automatic' decoupling in the RGE's; namely, it enjoys simultaneously 
that 

(i) the quadratic running of the scalar mass parameters is absent, 
(ii) the decoupling effects of heavy particles are naturally built in. 
Based on these properties (i) and (ii), we show, by adopting a simple model with 

[4] 

two mass-scales, that the same condition as in the single-mass-scale case is 
enough to achieve the RG improvement of the effective potential over the whole 
region of field space. 

We should remark that the property (i) is crucial to prove the above statement. 
Generally in MD scheme, there appear non-logarithmic and power-like corrections 
proportional to /i 2 , which are potentially large in the high-energy region. Such 
non- logarithmic corrections cause trouble in summing up the leading logs. We 
modify the renormalization scheme in order to cure this point. 

The existence of non-logarithmic corrections is related to the scheme depen- 
dence of the RG improved potential. [To examine this point is another motivation 
of the present work.] Note that it is not trivial at all that the RG improved poten- 
tials in MS and MD schemes coincide with each other. Of course, the full effective 
potential is independent of the renormalization scheme: the effective potentials in 
various schemes are related with each other simply by changes of variables. The 
effective potential correctly calculated up to a certain loop order is also scheme in- 
dependent since the loop expansion has a scheme-independent parameter, Planck 
constant h. In general, however, once one makes an approximation to the full 
theory, it is quite possible that the results are different scheme by scheme; some 
schemes give better approximations than the others. 
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In our case, we approximate the full effective potential by resumming 'loga- 
rithmic' parts of the perturbation series so that it satisfies the RGE. Then the 
scheme independence becomes nontrivial: the RGE relates 'log factors' at different 
loop orders, but the 'logarithmic' structure will differ scheme by scheme. More- 
over, there may appear non- logarithmic corrections as mentioned above. This is 
why there is no a priori relation between the RG improved potentials in various 
schemes. Do they give the same approximations? This is the problem of the scheme 
dependence. Our result will support to some extent the naive expectation that it 
is scheme independent. 

This paper is organized as follows. In the next section, we briefly review the 
basic ingredients for improving the effective potential by the RGE. In sections 3 
and 4, we define the MD renormalization scheme and discuss its basic features. 
We show that the RGE's in MD scheme inherit the nice property (i) as in mass- 
independent (MI) scheme ^ as well as the property (ii) as in MD one. The absence 
of the quadratic running is proved directly from the renormalization conditions and 
the automatic decoupling is established by utilizing the decoupling theorem. The 
detailed study on the structure of the effective potential in MD scheme is given 
in section 5. We first define the leading log series expansion in MD scheme and 
describe how to sum up the leading log. It will be shown, by examining the high- 
and low-energy regions separately, that we can correctly sum up all log factors 
over the whole region and that non- logarithmic corrections are in fact small. After 
establishing the procedure in MD scheme, we compare the leading log potential in 
MD scheme with those in other schemes such as MS and MD ones in section 6. By 
numerically solving the RGE's, a good coincidence will be found. A final section is 
devoted to conclusions and further comments. Some one-loop results can be found 
in appendices. 
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2. Improving the Effective Potential in MS scheme 
and Problems of Multi-Mass-Scales 

In order to explain the basic ingredients needed later, let us first make a review 

[4 7] 

of the procedure ' for improving the effective potential by using the RGE in MS 
scheme. We also describe why the problem of multi-mass-scales arises in the context 
of RG improvement of the effective potential. 

Following Ref. [7], let us consider the Yukawa model 

£ = \^? - \m 2 4> 2 - ^A0 4 + - g<f>)ip - u (2.1) 

where is a massive real scalar field and ip — (4>i, • • • ■,' 1 Pn) T are massless Dirac 
fields. We take the Dirac field to be iV-component in order to indicate which 
correction comes from fermion loop. For simplicity, we impose 'chiral-parity' in- 
variance, — > — 0, ip — > ^ip, to forbid the bare mass of fermion. The last term 
uj{= hm 4 in the notation in Ref. [4]) is a vacuum-energy term, which is usually 

[4] 

omitted but plays an important role in MS scheme. In this paper, we assume 
that both coupling constants g 2 /(167r 2 ) and A/(167r 2 ) are small and of the same 
order. 

In order to compute the effective potential V((p) for the scalar field VEV, 
<p = ((f)), we make a field shift — > + ip in (2.1), and obtain 

C =i(90) 2 - -M^ 2 - -M<f - ^A0 4 

2 V YJ 2 bY 3! ^ Y 4! Y (2.2) 

+ — M F — (70) ip — V^\(p) + [ 0-linear terms ] 

where the last term is the tree potential V^°\<p) = uj + m 2 ip 2 /2 + A(/? 4 /4! and we 
have introduced the masses for the boson and fermion ip, respectively, 

Ml =m 2 + \p 2 , M F =gip (2.3) 

in the presence of the scalar background ip. When ip is small, the field may be 
regarded as heavy field and ip as light field. 
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The starting point is that the effective potential is independent of the renor- 
malization point ji and thus satisfies the RGE 



= fi-^V = VV(<p, A, g 2 , m 2 , u ; , (2.4) 

d d d d d d 
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supplemented with the RGE's for parameters, such as 



[i^-g 2 = P g , /i-^-m 2 = -7 m m 2 , p-j-u = -Icum 4 . (2.6) 



One can immediately write down the general solution to (2.4) as 



V(<p, A, g 2 , m 2 , u;n) = V(?(t), X(t),g 2 (t),m 2 (t),uj(t) ; ety (2.7) 

where the barred quantities (p(t), etc., denote the solutions of running equations 
with a running distance t from the initial values ip, etc., at the renormalization 
point /j. [Here we are regarding the RGE's as differential equation with respect to 
an independent 'time' t, not /i] 

The RG improvement of the effective potential consists in solving the RGE 
(2.4). The RGE (2.4) by itself, however, does not determine the RG improved 
effective potential since it is the first order homogeneous differential equation. We 
should impose the suitable boundary condition on the functional form of V at a 
certain 'time' t. We call the boundary 'value' of the potential boundary function. 
The RG improved potential is fixed by requiring that the R.H.S. of Eq. (2.7) coin- 
cide at a certain 'time' t with the boundary function. [The RGE (2.4) guarantees 
that we can make a convenient choice of t.] It is the choice of boundary functions 
that determines how well the obtained potential approximates the exact one. 

How can we find a suitable boundary function ? 
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Let us work in MS scheme for a moment and examine the detailed structure 
of loop corrections to the effective potential. [The following arguments are valid 
in any mass-independent (MI) schemes since the structure of loop corrections does 
not change.] The simplest way' 4 ' to see this is to rewrite our Lagrangian (2.2), by 
rescaling the fields by a factor g as $ = g<p and \& = gip, into 



g 2 



-(d$) 2 - -M 2 ® 2 + - M F )^ 

2 2 , 

1 A 1 A i (2 ' 8) 



and to regard $ and \P as our basic quantum fields. In this form, the parameter 
g 2 is an overall factor in front of the action just like Planck constant K. So, L-loop 
contribution (L > 1) to the effective potential clearly takes the form: 



, -(/ ) 2 4 / 9 2 \L [, • , M 2 M 2 M 2 A i , n oN 

V = • (^) x [ function m In -i ln-f, ^, ? j . (2.9) 



In MS scheme, we have two types of logarithms ln(M 2 //i 2 ) and ln(M 2 //i 2 ) in 
our two mass-scale system, both of which can become large. Since we know the 
logarithms appear at most to L-th power at L-loop level, we can rewrite Eq. (2.9) 
into 

.,4 i+j<L 2 

^ (L) = §- E (^) L " (i+i) i L) (x,y)44 (2-10) 



tj>0 



by introducing the variables 

a 2 M 2 a 2 M 2 

Although we are assuming that the coupling constant g 2 / (167T 2 ) is small, we should 
regard the 'Kastening variables' s F and s B as 0(1) since the logarithms may be 
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large. Other variables x, y and hence the coefficient functions v\j (x, y) are also 
0(1). Then we sum up with respect to L and further rewrite it into the 

summation over £ = L — {i + j) as 

^=E^>=" + ^fE (i^)'A. (2.13) 

L=0 y ^=0 

oo 

/ £ ( SF , SB ;x, y ) = ^ +i+j) (^)44 ( 2 - 14 ) 

where we have included the tree part, fo = x^ 1 /2 — 5y/24, into the summation. 
This form of the expansion of the effective potential, first introduced by Kasten- 

[31 

ing for the single mass-scale case, is called leading log series expansion. When 
expressed in terms of the variables (2.11) and (2.12), it is the power series expan- 
sion in the small coupling constant g 2 /(IQn 2 ). The coefficients fo, fi, - • • , ft, - • • 
correspond to the leading, next-to-leading, • • •, £-th-to-leading, • • • log terms, re- 
spectively. Of course, it does not matter whether one uses A instead of g 2 as the 
expansion parameter. 

Now let us return to the question of how to specify the boundary function. 
The summation in (2.14) for the £-th-to-leading log term fi involves the quantities 
at L — £, t + 1, • • • loop level. If one could set s F = s B = in (2.14), then only the 
first term with i — j — would survive and the summation would ternimate at 
finite loop order, L — £, 

f e (0,0;x,y) = v { L =£ \x,y) . (2.15) 

Since Vqq can be obtained by computing L-loops, this would imply that one could 
use the £-loop potential, Vg = + • • • + evaluated at s F = s B = as 
the boundary function for the £-th-to-leading log potential. In other words, if one 
could find a 'time' to such that 

S P (i ) = S B (*o) = , (2.16) 
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then the desired £-th-to-leading log potential would be given by requiring that 
the R.H.S. of Eq. (2.7) be the boundary function Vl = i\ Sf=Sb= q: 



Actually, the condition (2.16) is sufficient, but not necessary one. For our 
purpose, it would be enough to find a 'time' to at which the logarithm factors 



As can be seen from Eq. (2.10), under this condition, these log factors contained at 
L-loop level reduce to precisely L-th-to-leading log order quantities. So, to obtain 
the boundary function for £-th-to-leading log potential, it would be necessary and 

[71 

sufficient to retain these log factors up to £-loop, just as in Eq. (2.17). 

Unfortunately, such condition (2.16), or even weaker one (2.18), can not always 
be satisfied simultaneously since the difference s B — s F becomes 0(1) when g 2 Lp 2 <C 
m 2 . The RGE (2.4) enables us to set just a single variable to desired value, but 
not several variables. So one can not find a solution to (2.16) or (2.18) and is left 
with the infinite summation in (2.14). This is the problem of 'multi-mass-scales' 
in the context of the RG improvement of the effective potential. 

Now, we examine whether MD scheme provides us with a solution to this prob- 
lem. Even in MD scheme, one will have a similar structure of the leading log series 
expansion as in Eq. (2.13). However, there explicitly appears the renormalization 

* As was proved in Ref. [4], one should use the RGE's at (£ + l)-loop order. Note also that, 
strictly speaking, the error in this equation is 0((g 2 /K>ir 2 ) i+1 ), not 0((g 2 /167r 2 ) £+1 ), but 
the difference will be small unless some coupling blows up (where our approximation itself 
does not make sense). 



V(<p, A, g 2 , m 2 , u;n) = V L=e ((p(t), X(t), g 2 (t),m 2 (t),u)(t) ; e%) 



S F (*)=0 

s B (t)=o 




(2.17) 



s F (to) and s B (to) are of 0(g 2 /16-7T 2 ), instead of zero, 




(2.18) 
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point [i in the effective potential, other than In Such explicit [i dependence is 
closely related to the existence of the quadratic running of scalar mass and makes 
the perturbation theory unreliable. In particular, in the context of the leading log 
series expansion, it may make the coefficient functions arbitrarily large. In 
fact, as we shall show in the following sections, we remedy this point by modifying 
the renormalization conditions in MD scheme. 



3. Modified Mass-Dependent Scheme 

In this section, we give a definition of the modified MD scheme (MD scheme). 
For the theory (2.1), we define the MD renormalization scheme by the following 

(2) 

renormalization conditions. For the scalar two-point vertex T\ , we impose 



.(2) 



p =0 
m 2 =0 



= lim 

m 2 ^0 



9 p(2) 



dm 2 



= -1 , 



0, 



9 r(2) 



dp 1 



(3.1) 
(3.2) 
(3.3) 



3 2 =— /i 2 



The fermion two-point vertex takes the form = A(p 2 )]f> due to the 'chiral- 
parity' symmetry, for which we require 



A 



p2 — _^2 



(3.4) 



Finally, for the Yukawa vertex and the scalar four-point vertex r^ 4 -*, we impose 



rf( P , -p;0) 



p 2 =— n 2 



-g, r 



(4) 



PiPj=-H 2 5ij+\li, 2 {l-5ij) 



-A (3.5) 



(3) 

where we set the boson-external momentum equal to zero in Tg . 
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To be precise, these renormalization conditions (3.1)- (3.5) should be supple- 
mented with that for the zero-point vertex 



r (0) = JL r (0) 
m 2 =o dm 2 



m 2 =0 



Clearly, the vacuum-energy term to = hm 4 , which played an important role in MS 
scheme, is completely independent of the renormalization point and is irrelevant 
for later discussions. [Instead, one can simply impose = —uj.\ 

The new set of the renormalization conditions (3.1)-(3.5) is a modified version 
of MD ones. The modifications are made in the conditions on the scalar two-point 
vertex, Eqs. (3.1) and (3.2), which take the place of a single condition 



r (2) 



fi 2 - m 2 (3.7) 



p2 — _^2 



in the conventional MD scheme. With this modification, MD scheme enjoys the 
properties announced in the introduction; to be precise, 

(i) In the high-energy region /i 2 ^> m 2 , the RGE's in MD scheme approach 
to those in a certain mass-independent (MI) scheme. In particular, the mass 
parameter runs at most logarithmically; 

(ii) In the low-energy region /i 2 <C m 2 , the decoupling effects are automatically 
taken into account in the RGE's and the vertex functions. 

As we shall see, it is crucial to separate the condition (3.7) into Eqs. (3.1) and (3.2) 
in order to realize the property (i), which will play important roles in section 5. 

Before showing the properties (i) and (ii) , let us take a close look at the renor- 
malization conditions (3.1)-(3.5). First, they differ from those in MI scheme! 1 ^ 
In MI scheme, one treats one parameter family of theories with different values of 
mass and renormalizes them at certain value like m 2 = 0, m 2 = /i 2 ; one imposes, 
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(2) 

for instance for 1^ , the condition (3.1) and 

9 -,(2) 



9 p(2) 



<9m 2 



- -1 — r v 

p 2 =V ' dp 2 * 

m 2 =0 



2 2=1- (3.8) 

p =-n 

m 2 =0 



Clearly all the renormalization constants are independent of the renormalized mass 
parameter m 2 . On the other hand, in our MD scheme, we are still treating, in a 
sense, one parameter family of theories with different values of mass in order to 
impose the conditions (3.1) and (3.2). With the renormalization conditions (3.2)- 
(3.5), however, the renormalization constants Z x for X — (f>, tp, m, g and A 
generally depend on the ratio fj?/m 2 ; 



Z x = Z x (\,g 2 M^;^) . (3.9) 

The RG coefficient functions {j3 and 7 functions), which are calculated from Z's, 
also depend on the mass parameter. 

Secondly, the renormalization constants Zx are consistently determined in MD 
scheme. A complication occurs only in the scalar two-point vertex while other 
vertices can be treated in the same manner as in MD scheme. Let us write the 
scalar two-point vertex as 

Vf(p, -p ; m 2 ) = Z^p 2 - Z m m 2 + Tl(p 2 ; m 2 ) . (3.10) 

As usual, the wave-function factor is determined by the condition (3.3): 

(3.H) 

As for Z m , the renormalization conditions (3.1) and (3.2) yield, respectively, 

lim m 2 Z m — n| 2 n , (3.12) 



1 9 n 

dp 2 



p 2 =—fj? 
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( 1+ - , p)«-= 1+ [( 1 -4"p n ) 



(3.13) 



p 2 =—y? 



where we have used Eq. (3.11) in deriving Eq. (3.13). Observe that Eq. (3.13), being 
a differential equation, does not completely determine Z m . This is most evident 
by noting that a piece Z m ~ /i 2 /m 2 drops from the L.H.S. of Eq. (3.13). What 
determines this piece is precisely Eq. (3.12). Thus we see that the renormalization 
condition (3.1) provides the condition (3.2) with a boundary condition and that 
all the renormalization constants are uniquely determined in MD scheme. 

Now, let us look at one-loop examples and confirm that the properties (?) and 
(ii) actually hold. [See appendices for more details.] The (3 and 7 functions are 
given by 

2 2 

16n\ = 2Ng 2 , 16n 2 ^ = > 

16n 2 lm = -ANg 2 , 16vr 2 ^ = 6g*K g (^) + Wg* , (3.14) 

2 

IGttV = , 16tt 2 /? a = 3AW^) + 8Ng 2 (X - Qg 2 ) . 

The functions Kx(z) (X — ip, g, A) are defined by 





2 


2 


1 + 








z 


z 


1 + 


2 


4 






Yz 


~ Yz 




3 


1 



1 

z 

1 



^ = 1 + ^-^(1 + ^)^ + 1), (3.15) 
K xi z) = l-- / ' ln V^+! 

which are normalized to be 1 in the high-energy limit z (= yU 2 /m 2 ) — > 00 and, 
remarkably, vanish in the low-energy limit z — > [See Fig. 1.]: 

f 1 as z — > 00 
X x (z) - . (3.16) 

1 as z^O 

Recall that the terms proportional to iV come solely from the light-particle 
(fermion) loops. Others come from the heavy-particle (boson) loops. The lat- 
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ter terms are accompanied by the functions Kx, which have the property (3.16). 
This is nothing but the decoupling of heavy particle loops, as claimed in (ii). 



Fig.l 



In MI (or MS) scheme, the RGE's do not have such a property of the automatic 
decoupling. Instead, one has to switch from the full theory to the low-energy 
effective theory. In MI scheme* the RGE's for the full theory are 

a 2 

16tt 2 70 = 2Ng 2 , I6n 2 ^ = y - , 

167r 2 7m = -ANg 2 , 16vr 2 ^ = 6«? 4 + ANg 4 , (3.17) 
167r 2 7„ = ~ , 16tt 2 /? a = 3A 2 + 8Ng 2 (\ - Qg 2 ) . 

In the low-energy effective theory, we keep only the terms proportional to N in 
Eqs. (3.17) and have the RGE's 

16tt 2 70 = 2Ng 2 , 16n\ = , 

167r 2 7m = -ANg 2 , 167T 2 = 4Ng 4 , (3.18) 
16ttV = , 16tt 2 /? a = 8iV£ 2 (A - Qg 2 ) . 

By comparing Eqs. (3.14) with Eqs. (3.17) and (3.18), one clearly sees that the 
RGE's in MD scheme interpolate those in MI scheme for the high- and low-energy 
regions. 

At one-loop order, the 7 m in MD scheme is the same as in MI scheme: 

d 2 4i\y 2 

^ m = ' (3 - 19) 

which means that the mass parameter in MD scheme runs logarithmically, as 



(2) 

* Here we adopt the renormalization condition as in Eq. (3.8). If we renormalize at 
m 2 = fi 2 , instead of m 2 = 0, then the -f m in such MI scheme coincides with that in MS 
scheme: 167r 2 7 m = —ANg 2 — A. But the difference is not so important here. 



-14- 



claimed in (i). [See Fig. 2.] This is the result of our modification of the renormal- 
ization conditions. 



Fig-2 

This is in sharp contrast to the case of the conventional MD scheme. Indeed, with 

(2) 

the MD renormalization condition (3.7), a fermion one-loop contribution to 
produces ^ 15 ' 16 ^ a piece proportional to /i 2 /m 2 in 7 m ; 



^ = 16^ (m (3 - 20) 

One sees that the running of the mass parameter is completely different from that 
in Eq. (3.17) since the second term on the R.H.S. of Eq. (3.20) dominates in the 
high-energy region (while it approaches to that in Eq. (3.18) in the low-energy 
region). If such quadratic running is present, the RGE's never interpolate the MI 
ones in the high- and low-energy region. 



4. Logarithmic Running and Automatic Decoupling in RGE 

In the last section, we illustrated the properties (i) and (ii) of MD scheme 
by one-loop examples. We now present the general argument to show that these 
properties hold to any loop order. 

Basically the property {%) follows from the fact that we introduce the \i de- 
pendence only through the dimensionless combinations of vertex functions r( n ) 
(n 7^ 0), such as (d/dm 2 )T^ 2 \ In other words, we never introduce the \i de- 

(2) 

pendence in the renormalization condition (3.1) on which has dimension two. 
Since the dependence on \i is introduced only through quantities which are at most 
logarithmically divergent, we do not meet the quadratic dependence on \i. 

This property {%) can be confirmed directly as follows. First, the condition 
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(3.3) determines the momentum dependence of the two-point vertex to be 



,(2) 



,m 



2\2 



2 2 2 /' \ , 



(P 2 + A * 2 ) 



9 9 

-p m 



(4.1) 



The condition (3.2) implies that the unknown function c is independent of m 2 , and 
the condition (3.1) determines it to be /(0, 0). Note that the function / should 
be nonsingular in the limit m 2 — > since has a massless limit. Then the 
renormalized scalar two-point vertex in MD scheme takes the form 



p + fj> 

~7 2 



2 x2 



-2 _2 



(4.2) 



Now we use the RGE for T 



(2) 



(£-2 7 ,- 7m ^K) 



dm 2 
d 



(4.3) 
(4.4) 



Inserting the general form (4.2) into Eq. (4.3) and taking a limit m 2 — > after 
setting p 2 = 0, we obtain that 



= lim m 7 m 

m 2 ^0 



2 9 _, n m 2 



(4.5) 



Since the quantity in the square bracket does not vanish (at least perturbatively) , 
Eq. (4.5) implies that 



= lim m 7 m . 

m 2 ^0 



(4.6) 



Since the quadratic running of the mass parameter corresponds to the behavior 
lm ~ fi 2 /m 2 , Eq. (4.6) proves the absence of the quadratic running in MD scheme. 
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Next, we turn to the property (ii). Let us examine the relation between the 
full theory in the low-energy region and the low-energy effective theory. Here the 
low-energy effective theory is obtained from the full theory by regarding heavy 
fields as external fields (instead of quantum fields), i.e., by taking out heavy-field 
internal lines. What we want to prove is that in MD scheme, the full theory in the 
low-energy region will automatically go over into the low-energy effective theory. 

In order to find such relation, we make use of the decoupling theorem: ' 8 ' 9 ' 10 ' 
the contributions due to heavy particles, aside from those which are suppressed by 
the inverse power of the heavy mass, can be renormalized into the parameters of 
the low-energy effective theory. Let be n-point vertex in the full theory and 
r( n ) the corresponding vertex in the low-energy effective theory. [We denote the 
quantities in the low-energy effective theory by the tilde.] Then, according to the 
decoupling theorem, when all the external momenta pi as well as the renormaliza- 
tion point fx are small compared with the mass m, i.e., for \piPj\, fJ? <C m 2 , 

T^( Pt ,g,X, m 2 ^) = ^4fW( K) P,m 2 ;,) + 0(^^) (4.7) 

where b and 2f stands for the number of external bosons and fermions, respectively: 
n = b + 2f. r( n ) does not depend on m 2 except for n = 2. 

Originally, the low-energy effective theory is not completely fixed by specifying 
the Lagrangian itself. So we fix it by imposing the same MD renormalization 
conditions as (3.1)-(3.5). Then, let us look at the scalar two-point vertex 

Tf(p,-p) = Z^f(p,-p) + 0{^-,^) x[-p 2 or» 2 } (4.8) 

where we have retained a factor — p 2 or ji 2 . [m 2 never appear here.] By differenti- 
ating this equation with respect to p 2 and setting p 2 = —/j 2 , the conditions (3.3) 
for T { ? and f { ? lead to 



l+0(£). (49) 
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Similarly, one can use the conditions (3.4) to show that 

^ = l + o(i£). (4.10) 
Thus, the relation (4.7) between r( n ) and r( n ) reduces simply to 

r<»> (pi, g 2 , A, m 2 ; „) = f W ( p^g 2 , A, m 2 ; „ ) + O (™>- 4 ) ■ ( 4 - n ) 

As for the dimensionless couplings g and A, we set p 2 = — /i 2 in Eq. (4.11) with 
n = 3, 4 and use the conditions (3.5) in the full and the low-energy effective theories 
to obtain 

9 =5+o(^), A = A + (^). (412) 

In this way, the finite renormalizations are not necessary also in the coupling con- 
stants. It remains to show that the same is true for the mass parameters: 

2 

m 2 = m 2 + 0(/i 2 ) = m 2 x l + o(J^j ■ ( 4 - 13 ) 

(2) 

We use the general form of the scalar two-point vertex 1^ , (4.2), from which we 
have, by setting p 2 = —/j?, 

= V-^ 2 -^ 2 /(0,0) . (4.14) 



r (2) 



~(2) 

We also have the same expression for '. Inserting both the expressions into 
Eq. (4.11) with n = 2 and p 2 = —fj?, we have 

2 

_ m 2 _ Ai 2 /(Q> Q) = _^2 _ ~ 2 _ ^2 J (Q) 0) + ^2 ^ 

which is nothing but the desired result (4.13). Thus, once the low-energy effective 
theory is renormalized by the same MD conditions, we no longer need the finite 
renormalization relating the full theory in the low-energy region to the low-energy 
effective theory; all parameters in MD scheme automatically go over into those in 
the low-energy effective theory. 
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Since the parameters in high- and low-energy theories are related in a way 
described above, it is easy to see that the RG coefficient functions in both theories 
are the same modulo 0(/j?/m 2 ) corrections: for X — g, X and Y = 4>,ip,m 2 , we 
have 

Px=(3 x + 0(^), 7r = 7r + 0(^). (4.15) 

This completes the proof of the property (ii). 

It is instructive here to see how the MD scheme modifies the conventional 
MD one. Let us apply the same reasoning as above in the proof of (i), to the 
conventional MD scheme. Again, due to the condition (3.3), the renormalized 
two-point vertex takes the form (4.1). Now, the renormalization condition (3.7) 
determines the unknown function c(g 2 , A; m 2 //i 2 ) to be zero: 

Tf MD = P 2 -m 2 + 0((p 2 + ^)- (4.16) 
Inserting this into the RGE (4.3) and setting p 2 = — /i 2 , we obtain 

2 

7m = -2(l + £2)7*. (4-17) 

This clearly shows that the mass parameter in MD scheme runs quadratically in 
the high-energy region (as long as the wave-function renormalization 7^ does not 
vanish) . 

Let us recapitulate this in some different way. In MD scheme, the scalar 
two-point vertex (4.16) satisfies all the MD renormalization conditions except the 
condition (3.1), and we meet the quadratic running of the mass parameter m 2 D . 
Now, we finitely renormalize the two-point vertex so that the condition (3.1) is 
satisfied. By requiring the equality between the expressions (4.16) and (4.2) at 
p 2 = —fi 2 , we have 

™L = ™ls + A* 2 /(0,0) . (4.18) 

This shows that the mass parameter m 2 ^ is just the logarithmic part of m 2 D : 
we have succeeded in separating the logarithmic and quadratic parts in m 2 D by 
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dividing the MD renormalization condition (3.7) into the MD ones (3.1) and (3.2). 
This is the property (i). 

In this way, we conclude that the present MD scheme simultaneously enjoys 
the 'automatic' decoupling of heavy particles (as in the conventional MD scheme) 
and logarithmic RG running (as in MI scheme). 



5. Improving the Effective Potential in MD scheme 

We now turn to our main task of how to improve the effective potential in 
the presence of several mass scales by using the RGE in MD scheme. Let us first 
describe the structure of the effective potential for the system (2.8) in MD scheme. 
This can be done by applying almost the same reasoning as reviewed in section 2 for 
MS case. A difference arises from the finite part of counter-terms:* there appear 
another type of log factor ln(m 2 /// 2 ) and non-logarithmic dependence on /i. So 
L-loop contribution now takes the form 



* r( j,\ A/* , g 2 \L[ r , , . M 2 Mg m 2 

V y ' = — 7y-{ 7y ) L-tn order polynomial m in — In — in — ~- 

gZ ^\Q n z; [ n z n z 

whose coefficients depend on 



M 2 ' g 2 ' m 2 



(5.1) 



Introducing the variables 



a --^-In^ s --^-In^ s =^ln— (5 2) 

SF ~ 167T 2 ^ /i 2 ' SB " 167T 2 ^ ^ 2 ' 167T 2 M /X 2 ' (5 ' 2j 

X = M' y = ?' z = (5 - 3) 



* Loop corrections themselves depend on m 2 only through the combination M 2 = m 2 +\ip 2 /2 
since we evaluate the potential in the background ip = ((f>). But since the renormalization 
constants are determined in the symmetric phase ip = 0, counter-terms produce the depen- 
dence on the ln(m 2 /^ 2 ) and [i 2 /m 2 . 
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we rewrite Eq. (5.1) (L > 1) into 

i\/r4 i+j+k<L 2 

^-f £ (i^) L - (i+j+k) 4%(z.y,z)44si- (5.4) 

i,j,k>0 

Then we define the leading log series expansion in MD scheme by 

f£(s F ,s B ,s m ;x,y,z) = ^ vf^ +J+k \x,y, z) sls{s k m (5.6) 

i,j,fe>0 

where we have again included the tree part into the summation and we have defined 
the order t of this expansion as t = L — (i + j + k). 

The leading log series expansion (5.5) is the power series expansion in a small 
coupling constant g 2 /(167r 2 ), as before. Compared with (2.14) in MS scheme, 
however, the coefficient functions have the dependence on the new variable 
z = /i 2 /m 2 which potentially makes large. So we do not know, at this stage, 
whether or not the expansion (5.5) is sensible one; we do not know whether or not 
terms in the (I + l)-th-to-leading log order are smaller than those in the £-th-to- 
leading log order. Furthermore, we have a log factor s m ~ ln(m 2 / /i 2 ) in addition 
to the 'original' ones (2.11), s F and s B . So it is not obvious how one can sum 
up these logarithms simultaneously by using the RGE which can eliminate just a 
single variable. 

At first sight, the situation in MD scheme appears to be worse than in MS 
scheme due to the new variables s m and z. What save the day are the nice prop- 
erties established in the last section. Based on the properties (i) and (ii) as well 
as the fact that the model has the well-defined massless limit, we claim that the 
single condition 

S P (t) = (5.7) 

is enough to determine the correct boundary function over the whole region of field 
space. The potentially dangerous variable z = /i 2 /m 2 are in fact harmless and the 
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remaining logs s B and s m are automatically summed up; otherwise they decouple? 
As a result, with £-loop effective potential Vl = i and (£ + l)-loop RGE's in MD 
scheme at hand, the £-th-to-leading log potential is given simply by 

V(<p,X,g 2 ,m 2 ,u;-fi)^V L= i(<p(t),X(t),g 2 (t),m 2 (t),o(t)-e t pL) 

Sf (t)=0 

t ^«of(^n. (5 ' 8) 

The proof of the statement proceeds in a regionwise manner; we divide the field 
space into the large <p (high-energy) region g 2 (p 2 ^> m 2 , the small (p (low-energy) 
region g 2 (p 2 <C m 2 and the intermediate region g 2 Lp 2 ~ m 2 , and prove that the 
statement holds region by region. We should stress here that the final answer (5.8) 
does not require to divide the region of field space unlike the procedure in MS 
scheme. 

First, the proof is rather simple in the intermediate region g 2 (f 2 ~ m 2 . In this 
region, s B and s m are of 0(g 2 /I6ir 2 ) from Eq. (5.7) and the weak condition as in 
Eq. (2.18) is satisfied. So, s B and s m are already summed together with s F by the 
single condition (5.7) and the correct choice of the boundary function is the same 
as in the single mass-scale case, i.e., V^ = ^| sf=q . 

Concerning the asymptotic regions (g 2 ^ 2 ^> m 2 and g 2 (p 2 <C m 2 ), we make 
an observation needed for the proof. The boundary function is defined by setting 
/j 2 = M 2 (= g 2 i~p 2 ) and depends on m 2 only through the ratio m 2 /jj?. It follows 
that as far as the boundary function is concerned, taking the high-energy limit 
ji 2 = g 2 ip 2 — > oo is equivalent to considering the massless limit m 2 — > while the 
low-energy limit fi 2 = g 2 (p 2 — > corresponds to the limit m 2 — > oo. 

It should be noticed that the validity of this equivalence between the high- 
energy and massless limits heavily depend on the property (i). If m 2 run quadrat- 
ically, then the ratio m 2 / fj? would approach to a finite value in the high-energy 
limit, and the equivalence would break down. 



| This is why we have used the log factor of the lightest particle in the condition (5.7). 
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5.1. High-energy region (g 2 y? 2 > m2 ) 



Now, let us begin with the large cp region, g 2 ip 2 3> m 2 . The second log factor 
S B becomes asymptotically equal to the first logarithm S F 

*-*=i^[ to ^ +to ( i+ ^>]= o (isO (5 - 9) 

as is the case in MS scheme. [Recall y = X/g 2 = 0(1).] So, setting s F equal to 
zero is equivalent to setting s B equal to zero, modulo the quantities of 0(g 2 /I6ir 2 ) 
which is of the higher order in the leading log series expansion. Physically, we 
can regard the massive particle (here 0) as massless. Thus the condition (5.7) 
automatically sums up the second log factor s B as well as the first one s F . 

Next, let us make sure that the variable z = /i 2 /m 2 does not make the co- 
efficient functions large in the high-energy limit z — > oo. To this end, it is 
enough to see that there is no positive power term in z when the potential is ex- 
panded asymptotically in z~ 1 (<^i 1). As noted above, this limit is equivalent to 
the massless limit m 2 — > in our boundary function and we know that there arises 
no singularity in the latter limit. This establishes that power-like pieces of /i 2 /m 2 
vanish in this region. 

As for the third log factor s m , we show that it does not contribute to the 
boundary function. Since we know that the massless limit is regular, the dangerous 
variable s m in this limit disappears in our boundary function like 

9 9 

m . m 

—~ m — ~- — > . 

Thus, we no longer need to sum up this log factor in the high-energy limit g 2 (p 2 = 
/i 2 3> m 2 . This establishes our claim in the large (f region. 
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5.2. Low-energy region (g 2 (f 2 < m 2 ) 



Next, we turn to the small (f region. Now, s P = no longer implies s B = 
since their difference S B — S P becomes 0(1). Furthermore, s m is also large. So we 
can not sum up these two log factors, S B and s m , simultaneously. Instead, we shall 
show that in the low-energy limit g 2 (p 2 — > 0, the log factors s B and s m as well as 
z = fi 2 /m 2 decouple so that we no longer need to sum them up. 

This can be seen by using the decoupling theorem. As in Eqs. (4.12), (4.13) and 
(4.15), all the parameters as well as (3 and 7 functions of the full theory approach, 
in the low-energy limit, to those of the low-energy effective theory. In particular, 

2 

s, = s F + 0(^) • (5.10) 

This property of the automatic decoupling holds also for the effective potential 
itself. In particular, the boundary function l / | Sp=0 satisfies 

Here V is the effective potential in the low-energy effective theory fixed by the 
same MD renormalization conditions and we have replaced s F = with s P = by 
using the Eq. (5.10). [One-loop example can be found in appendix A.] From this 
expression, we see that there is no contribution from potentially large variables 
s B and s m as well as negative power terms in jj? /m 2 . The first term V does 
not contain loop effects due to heavy particle by definition. Furthermore, the 
remaining terms are small by itself. Thus, we need not sum up s B and s m . This 
establishes the correctness of our boundary function in the low-energy asymptotic 
region g 2 (p 2 = fi 2 <^ m 2 . Together with the automatic decoupling (ii) in (3 and 7 
functions, this completes the proof of our assertion in this region. 

In this way, we establish the procedure to improve the effective potential. The 
final answer is given by Eq. (5.8). 
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6. Comparison with Other Schemes: Numerical Study 



In this section, we work with the leading-log order and demonstrate the results 
of the RG improvement. At the leading log order, we use the tree potential 
as the boundary function, and one-loop RGE's. 

Note that our procedure is such that, at individual point ip in field space, the 
value V((p) of the improved potential is evaluated by Eq. (5.8). Actually, this is 
not economical since we solve the running equations at each ip. As explained in 
Ref. [4], one can avoid this duplication by finding the value of <p corresponding to 
the solution to s F (t) = for each value of the running distance t. Namely, as we 
solve the running equations, we simultaneously obtain the value of the effective 
potential at 



2 2* 2 1 rv^o) 12 



(6.1) 



The last factor (p(0)/(p(t) does not depend on the initial value (p(0). We set fj, = m 
for later convenience. Then, Eq. (6.1) becomes 



g\t) p2(t) 

so that the region g 2 (p 2 ^m 2 corresponds to t^O. 



2t (6.2) 



We compare the improved potential in MD scheme with those in MS and 
MD schemes. For this purpose, we should match the renormalized parameters in 
different schemes in order to guarantee that we are treating the same system. For 
the leading log potential, the parameter relations should be exact also in the leading 
log order. By using the fact that the improved potential exactly satisfies the RGE 
(2.4) with one-loop (3 and 7 functions, we match the parameters at /i — m. Then, 
the parameter relations reduce to the tree-level ones since s m ~ ln(m 2 //i 2 ) = 0. 

We now show the result of our numerical calculations. We present it in a 
moderate case (a): g 2 = 0.55, A = 2.3 and in extreme cases (b): g 2 = 0.5, A = 2.5 
in which A blows up, and (c): g 2 = 0.7, A = 2.0 in which the vacuum becomes 
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unstable. We set m — 1 (as mass unit), N — 1 and the vacuum energy to uo = 
so that V(ip = 0) = 0. Fig. 3 shows the asymptotic behavior of the leading 
log potentials in MD, MS and MD schemes for the case (a). [Other cases are 
similar, so omitted.] The horizontal and vertical axes are ln(<^ 2 /m 2 ) and In |V/m 4 |, 
respectively. We find that the scheme dependence is quite small* This can be 
understood from the fact that since we are using the tree potential as the boundary 
function, the asymptotic behaviors in the large and small ip regions are mainly 
determined by the quartic coupling A and mass m 2 , respectively. 



Fig.3 



The behaviors of g 2 (t) and X(t) in MD scheme are shown in Fig. 4. The horizontal 
axis is the same as in Fig. 3. The scheme dependence is mild even in the extreme 
case (6) where the quartic coupling A blows up in the high-energy region and the 
case (c) where A becomes negative and the vacuum instability occurs. 



Fig.4 



Fig. 5 



In order to examine how much the RG improvement is obtained, we evaluate 
the difference of the improved potential V from the tree one normalized by 



(6.3) 



V(°)(<p) 

Figs. 5 show the results of Xmd, Xms an d Xmd, respectively, for each case of the 
parameter choices. As expected, the longer distance we run from t — 0, the larger 



* A rather good coincidence is found between MS and MD schemes. It is not clear to us that 
this persists to higher leading log order since the reasoning presented in section 5 will not 
apply to the conventional MD scheme. 
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improvement is obtained. The large improvement \x\ ~ 1 is obtained for 2t ~ 30 
since the variables s F and s B becomes 0(1) for our parameter choices. In extreme 
cases, even larger improvement is obtained, but the leading log approximation itself 
breaks down for these cases. 

Finally, we add a remark. We can further improve the approximation by re- 
quiring that the potential be correct not only in the leading log order, but also 

[4 71 

in one-loop level. ' For this combined approximation, we should use the one-loop 
potential as the boundary function while one-loop RGE's are enough. Also, the 
parameter should be matched at one-loop level by performing the finite renormal- 
ization at [i = m* 

7. Conclusions and Discussions 

We have discussed the issues concerning the RG improvement of the effective 
potential in the presence of several mass scales. Originally the coexistence of 
multi-mass-scales causes trouble in determining the boundary function needed for 
the general solution of the RGE. By adopting a simple model possessing two mass 
scales, we have seen in this paper that MD scheme provides us with a suitable 
choice of boundary functions without dividing the field space: the correct boundary 
function for £-th-to-leading log potential is just £-loop potential evaluated at s F = 0. 

The MD scheme is the new renormalization scheme proposed in this paper. Its 
crucial properties are the automatic decoupling of heavy particles and the absence 
of quadratic running of scalar mass. These properties enable us to show that the 
leading log series expansion is well-defined even in the presence of non-logarithmic 
corrections. Then, the procedure can be stated by a single condition over the entire 
region, which is the same as in the single-mass-scale case. 

We make a comment on other possible methods to handle the multi-mass-scale 

[7] 

systems. First, we already have the procedure in MS scheme. Such regionwise 

* By matching at fi = m, higher order terms are small in MS and MD schemes. This might 
not be the case in the conventional MD scheme. 
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procedure will be cumbersome especially when there are many mass thresholds. In 

MD scheme, calculations of the RGE's become harder, but once they are calculated, 

then various threshold effects are automatically taken care of. Our method will be 

more useful when some intensive investigations will be needed, such as scanning 

[171 

for large parameter space. On the other hand, Einhorn and Jones proposed 
to introduce several renormalization points \i{ by which the multi-log factors are 
simultaneously summed up. Their method is interesting, but the RGE's become 
partial differential equations. Our method described here involves solving ordinary 
differential equations, which will be much easier task. 

Finally we comment on possible applications of the method in this paper. 
Amongst, it is interesting to apply our procedure to the analysis of the Higgs 
potential in the supersymmetric standard model, in which many mass scales are 
present. When the supersymmetry breaking scale is rather high, we expect large 
improvement to usual analyses which make use of at most one-loop potential! 18 ' 19 ^ 
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APPENDIX A 

We gather one-loop results by the dimensional regularization. We first deter- 
mine the renormalization constants Zx = 1 + hZjP + 0(h 2 ) from MD renormal- 
ization conditions (3.1)-(3.5). We then calculate one-loop RGE's and the one-loop 
contribution to the effective potential. Our conventions are: d = 4 — 2e, 
1/e = 1/e — 7 + ln47r and Trl = 4. A mass scale jiQ is introduced so that the 
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dimensionality becomes correct. /jLq always appears as /ijj 6 J d d k/ [2^) d i and is iden- 
tified with the renormalization point fx in MS scheme. Otherwise, jiQ has nothing 
to do with \x which is introduced through the renormalization conditions. 

As usual, we renormalize the theory (2.1) in the symmetric phase, (0) = 0. 
First, the vacuum-energy at one-loop is 



l l 



4 167T 2 



m 



r 1 3 , m l 
- + - - In — 

e 2 ^ 



(A.l) 



From Wbare = ^ + hou^ + 0(h 2 ), the simplified condition r(°) = — uj leads to 



(i) 1 1 



4167T 2 



m 



1 3 

- + - + In 



(A.2) 



Second, the boson self-energy, T^?\p, —p) = Z^p 2 — Z m m 2 + H(p 2 ), is given by 



2, 2Ng 2 



2 

2 P~ 



i + 2 + ln ^ 



A/2 , 



e m z 



(A.3) 



The renormalization constant is determined as usual by Eq. (3.11) while the 
Z m is determined by solving the differential equation (3.13) and imposing the 
boundary condition (3.12) to be 



Z?> = 



2Ng_ 
16tt 2 



2 r 



1 



= + l + ln4 



(A.4) 



The fermion self-energy, Y^\p, —p) = Z^j> — E(p), is similar: 



7 (1) 



16tt 2 



l + ln4-^(4) 



m 



(A.5) 



where I${z) = (1 + l/z) 2 \n(z + 1) — (1/z) — 2. The vertex correction to the Yukawa 
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(3) 

coupling is, T y g (p, -p;0) = -Z g g + A s (p, -p ; 0), 



Ai 3) (p,-p;0) = ^ 



+ In ^-^Hf) 



m 2 



(A.6) 



where = (1 + 1/z) ki(z + 1). From the condition (3.5), we obtain 



,(i) 



16tH 



(A.7) 



The vertex correction to the quartic scalar coupling, T y ^' = —Z\\ + A\, is 
a little bit complicated. We separate the boson- and fermion-loop contribution, 
= A B + A F . Let pi be incoming external momenta (i — 1 - 4) and s = (p\ +P2) 2 , 
t = (pi + Pa) 2 and u = (p\ + P3) 2 . The boson-contribution consists of the s-, t- 
and -u-channel ones as A B (pj) = A B (— s) + A B (— t) + A B (— u) where 

i2/o r 1 ,,2 / _q\ n 

(A.8) 



A.<-.) .^[I + ln*-7 A (^)l, 



= / da In 



1 H g a (l — 



In 



— s 
4m 2 



4m 2 



+ 1 — L 



Am? 
1 + 1 



(A.9) 



where L(£) = £ (hiC — !)• The fermion-contribution is evaluated in appendix B: 
Ng 4 



Jo(Pi) = 4 + 2 + In 



Jo(Pi) + [ permutation in pj ] , (A. 10) 

2 + J(-pf, -p| ; -(P1+P2) 2 ) 



t4 



-(pi +P2) 



\ J ( p\p\ , p\p\ ; (pi + k) 2 (p2 + ps) 2 ) 



(A.11) 



and the function J(£,r)] () is defined in the appendix B. Both A B and A F are 
completely symmetric in p\ {% = 1-4). Then, we impose the condition (3.5) at the 
symmetric point, pf = — /i 2 and s = t = u = — (4/3)yU 2 , to obtain 



A lg7r 2 



+ 1u4-/a( 



3m^ 



4!A^ 
16tt 2 



^ + ln4 + ^ 



(A.12) 



with a constant F = 3.02198 • • •. 
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(1) 



Some remarks are in order. First observe that functions Ix (contained in Z x 



for X — ip,g, A) behave in the high-energy (massless) limit /i 2 /m 2 — > oo as 



1*1*, 



In 



if! 



2 d (n 



in 



1 . 



(A.13) 



An important point is that in the low-energy limit /i 2 /m 2 — > 0, they all approach 
to constant values, 1^(0) = —1/2, I g (0) = 1 and I\(0) = 0, as 



^(5)-«0), „»<, x (£)^0. (A.14) 



Note also that we have to renormalize at the symmetric point in order that 
all the s-, t- and w-channel boson-loops equally contribute to I\. 

Now, we turn to the RG coefficient functions and the effective potential. The (3 
and 7 functions are obtained by noting that the bare parameters are independent 
of For instance, from g 2 w% = Z 2 Z^Z~ 2 g 2 , 



9 2 



ANg 2 



+ 



djj? 

2g 2 , d - 1 



24 1) +4 1) +< ) 



+ o(h 2 ) 



16n 2 16n 2 dji 



By looking at the asymptotic behavior (A.13) and (A.14), we introduce the 'inter- 
polating' functions (3.15) with the property (3.16) by 



d d Az 



2I g (z) + I^z) 



(A.15) 



With these definitions, we finally obtain the RGE's (3.14). 
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Next, we discuss the effective potential V(tp) for tp = ((f)). The one-loop con- 

.(i) 

loop 



tribution to the effective potential takes the form — vP^ + AV^ where 



V 



(i) 



loop 



+ 



64vr 2 



M 



1 


3 


+ ln 




€ 


~ 2 




1 


3 


+ In 




~€ 


~ 2 





(A.16) 



while the contribution from counter-terms are, from Eqs. (A. 2), (A. 4) and (A. 12), 



AF« ee a;* 1 * + zfi^mV + 4^ V 



+ 
+ 



4iV 
64tt 2 


M F 4 


fi+u.4 

L e /i z 


I 4iV 
J 64tt 2 




1 

64tt 2 


Mi 


- + '4 


+ ln-^ 
m z - 




1 

64tt 2 


3 4 i 2 \ 2 

—m + m \(f 
2 





(A. 17) 



Thus, the final form of V^ 1 ) in MD scheme is 
VKI) = 



64*2 r 



+ 



+ 



64vr 2 



m z 



3 

2 ~~ 2 

3 



1 



1 (-yf 



64tt 2 v 2 



In 



+ 



1 2 \ 2 



In 



V 2 



3i 
2 



64vr 2 
/4/i 2 



ln^ + 1 



AN 
64^ 



<? V ^ • 



(A.18) 



Note that the first two terms in Eq. (A.18) take just the same form as V^> in MS 
scheme while the remaining terms are extra contributions to it in MD scheme. 

From the expression (A.18), we can explicitly confirm what we have gener- 
ally argued in section 5. The expression (A.18) is suitable in the large f region. 
In the low-energy region /i 2 <C m 2 , we rewrite it, by combining ln(/i 2 /m 2 ) with 
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ln(M 2 //i 2 ), into the form 
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(A.19) 



where we used I\(0) = 0. When the second term is expanded in the small ip region 
(p 2 <C m 2 , terms proportional to m 4 and m 2 (p 2 cancel so that we have 



l ' ,1| = f ' ,1| + ^(V) 2 x°(4.^) 



(A.20) 



where we denote the first term in Eq. (A.19) as V^ l \ which is just the potential of 
the low-energy effective theory. This is an example of the 'automatic' decoupling. 

APPENDIX B 



,(4) 



We evaluate the fermion one-loop contribution A P to the four-point vertex 
— Z\\ + A B + A F . Generic one-loop integral can be reduced to the scalar- 
loop integral^ 20 ' and expressed in terms of Spence function, but we give another 
expression for A F * 



The fermion one-loop contribution to the four-point vertex is given by 



A P (pi,P2,P3,P4) = 



Ng* 



16tt 2 

I(P1,P2,P3,P4) = 167T 2 /io e 



I(Pi,P2,P3,P4) + [ permutation in p 2 , p 3 , p 4 ] , 



d d k 



Tr 



(B.l) 

where p — pi, q — pi + P2 and r = pi + p 2 + Ps = — Pa- This integral can be 
reduced to the scalar-loop integral by carrying out the trace and making use of the 
identity 2(k + a)(k + b) = (k + a) 2 + (k + b) 2 — (a — b) 2 for any momenta a and b. 
As a result, the integral / decomposes into three parts, I = I 2 + ^3 + 1 a, where I n 
contains n propagators. We denote as D a = (k + a) 2 . 



* The authors are grateful to T. Kugo for discussions on this calculation. 
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The first integral I2 is evaluated in cyclically symmetric way by making a 
suitable shift in the loop momentum as 



1 



d d k 


I 1 


(2vr) d i 


[ D p D r 


d d k 
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{2n) d i D D P1+P2 



+ 



1 



i4 



— + 2 + In - 

e -{P1+P2Y 



D D q 
+ [ cyclic in p { } 

+ [ cyclic in pi ] . 



(B.2) 



Similarly, the second one I3 is evaluated symmetrically as 



h = 16T J W)H D^D (> D K + 1 CyC " C 

= -P2 5 ~(Pi +P2) 2 ) + [ cyclic inpj ] . 



(B.3) 



We have introduced the function J(£, i]; () = J(rj, £ ; £) defined for £, 77, ( > by 



J(tV,C)=f da- ■ — + \ — 7^ N 

Jo £ a + r] [1 — a) — (, a{l — a) Qa[\ — a 



\J a+ ? (1 - a) (B.4) 



in terms of which we have useful formula 



., 1 d A k -a 2 -b 2 + (a + b) 2 , 2 2 2 , 

16tT / - — — = J( - a 2 , -6 2 ; -(a + 6) ) . (B.5) 



(27r) 4 i 
Finally, the ^4 is 



D- a D D b 



2 /" <i 4 /c r 2 (jo — g) 2 + p 2 (g — r) 2 — g 2 (r — p) 



2 „2/ 



(2vr) 4 i 



D DpD q D r 



(B.6) 



which can also be expressed in terms of J by making a conformal change^ 20 ' of 
the integration variable k^ = (/Iq/A; 2 ) k^. \jiq is an arbitrary scale.] Under the 
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conformal transformation = (/ig/a 2 ) a M , (A; + a) 2 = (a 2 /k 2 )(k + a) 2 , we have 



/• d 4 fc 1_ 

,/ [2;; s l; "71777"/) n ~~ n-',,-',^ 



r _d*h_ 
J I^) 1 



r) 4 i DoD p D q D r p 2 q 2 r 2 J (27r) 4 i DpDqDf 
By applying the formula (B.5), we evaluate the integral (B.6) as 
2 /" <i 4 £; — (p — §) 2 — (§ — f) 2 + (f — p) 2 



J 4 = 16tt z 



(27r) 4 i 



V g z r z r z £r / 

= --j- 7( PiPl , v\v\ i (pi +P2) 2 (P2 +Ps) 2 ) + [ cyclic in pi ] . (B.7) 

where we have used the definition (B.4) and the definition of p, q and r in the last 
equality so that the cyclic symmetry becomes manifest. 

Thus, the final form of the integral I = h + h + h and the fermion one-loop 
contribution A F are, from Eqs. (B.2), (B.3) and (B.7), 



I(pi) = 7o(pi,P2,P3,P4) + [ cyclic in pj ] , 



16tt 2 



7o(pi,P2,P3,P4) + [ permutation in pi 



(B.8) 
(B.9) 



with Iq given in Eq. (A. 11). At symmetric point, Jo reduces to the second term in 
Eq. (A. 12) with 



F = 2 + In^ + J(l, l;y)= 3.02198 



(B.10) 



where we numerically integrate the function J(£,r]; () = (£ + rj — C)Jo(£,,ViC)i 
(B.4), by using E(a) = £a + 77(1 — a) — C a (l — a ) an d 



+ da 
Jo E{a 



— j da 
'0 



1 £ a + rj (1 — a) 
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E(a) E(0) 
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In a + 



E(a) E(l) 



ln(l-a) 



(B.H) 
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FIGURE CAPTIONS 



Fig. 1 The 'interpolating' functions (solid line), K g (dashed line) and K\ 
(dot-dashed line). 

Fig. 2 A typical behaviur of m 2 in MD (solid line), MS (dashed line) and MD 
schemes (dot-dashed line). All schemes as well as the full and effective the- 
ories in MS scheme are matched at ln(/i/m) = 0. 

Fig. 3 The asymptotic behavior of the leading log potential in MD scheme (solid 
line), MS scheme (dotted line) and MD scheme (dot-dashed line). The inputs 
at [i — m are g 2 = 0.55, A = 2.3 (case (a)). V ~ ip 4 in the large (p region 
and V ~ (p 2 in the small ip region. 

Fig. 4 (a): The behavior of A (solid line) and g 2 (dashed line) in MD scheme. The 
input is moderate one: g 2 = 0.55, A = 2.3. 

(b) : Same as Fig. 4 (a) for the extreme case: g 2 = 0.5, A = 2.5. A hits 
Landau singularity in the high-energy region. 

(c) : Same as Fig. 4 (a) for the extreme case: g 2 = 0.7, A = 2.0. A becomes 
negative to cause the vacuum instability. 

Fig. 5 (a): The 'character' Xmb (solid line), Xms (dashed line) and Xmd (dot-dashed 
line) for the moderate case (a). The large improvement is obtained in asymp- 
totic regions. 

(b) : Same as Fig. 5 (a) for the extreme case (b). Due to K\, sizable scheme 
dependence is observed near Landau singularity. 

(c) : Same as Fig. 5 (a) for the extreme case (c). Since the mass term domi- 
nates the quartic term in V, sizable difference from MS and MD is observed 
in MD scheme. 
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Fig.3 
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